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ABSTRACT 



Following the formalism of Gell-Mann and Hartle, phenomenological equations 
of motion are derived from the decoherence functional of quantum mechanics, using 
a path-integral description. This is done explicitly for the case of a system interact- 
ing with a "bath" of harmonic oscillators whose individual motions are neglected. 
The results are compared to the equations derived from the purely classical theory. 
The case of linear interactions is treated exactly, and nonlinear interactions are also 
compared, using classical and quantum perturbation theory. 
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I. INTRODUCTION 



A. Decoherence and the quasiclassical limit 

Two of the most puzzling aspects of the quantum theory have, until recently, remained 
unclear: the proper interpretation of quantum probabilities, and the mechanism by which 
deterministic classical "laws" can arise from a probabilistic underlying theory. The idea of 
wave-function collapse, while providing a useful approximate description of most experimen- 
tal situations, begs the question of why a system which otherwise undergoes purely unitary 
evolution should suddenly and dramatically be collapsed upon measurement by a scientist. 
The procedure is highly asymmetric, instantaneous, and irreversible, and moreover requires 
the existence of a "classical" measuring device outside the system being measured. When 
one considers a closed system the idea of wave-functions collapsing becomes highly ambigu- 
ous. There is nothing outside the system to collapse it. The quintessential example of this, 
of course, is the universe itself. Clearly, if the fundamental laws of the universe are quantum 
mechanical, there can be no separate "classical domain" to explain our observations. Since 
the classical realm is itself, presumably, merely a limit of the underlying quantum reality, 
the probabilities must arise directly from the quantum theory itself, without recourse to the 
deus ex machina of the measurement device. And somehow, the various potential futures of 
the universe must collapse themselves onto the one possibility which we observe. Quantum 
cosmology requires a solid formalism for the treatment of closed systems, and work in this 
field should have that as its goal. 

The recent work on the decoherence functional formulation avoids the problems of earlier 
approaches []]-§]]. Physics is described in terms of exhaustive sets of possible histories, coarse- 
grained, with the restriction that these histories must be decoherent. That is, it must be 
possible to assign probabilities to these histories such that they obey the classical laws of 
probability, with no interference. 

Gell-Mann and Hartle have argued that it is possible, in a highly coarse-grained system, 
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to define the classical equation of motion directly from the decoherence functional itself 
0. I will, in this paper, attempt to show that this definition gives the exact classical 
results, at least for the case of systems interacting with baths of oscillators; and further, 
that these systems are decoherent in the classical limit. Quantum effects enter as a random, 
fluctuating force from the effects of neglected degrees of freedom, even in cases where the 
classical noise would ordinarily be zero. Fluctuations, dissipation, and decoherence turn out 
to be intimately interlinked. 

The linear case has been treated before by a number of people, both classically and 
quantum mechanically, though not in precisely this same way |6|-|8|]. The correspondence 
of this quantum system to the classical Langevin equation is thus nothing new. The deco- 
herence of similar systems has also been examined, using a somewhat different definition of 
decoherence which for these models generally corresponds to my definition ||. However, to 
my knowledge, no one has considered the classical correspondence of these sorts of nonlin- 
ear systems, nor the relationship between dissipation, noise, and decoherence in these more 
general cases. Thus, the results herein are of interest in demonstrating that it is possible to 
define classical equations of motion directly from the quantum theory in a broad range of 
systems. 

B. Path-integral description of the decoherence functional 

We will not, for the purposes of this model, be using the decoherence functional in its 
most general form. Instead, we will consider only one type of history. Suppose that our 
system is completely described by a set of generalized coordinates q 13 (collectively referred to 
as q). The most fine-grained possible family of histories would be just the set of all possible 
paths q(t). We can coarse-grain this by dividing the range of the q@ into an exhaustive set 
of intervals A 1 (U) at a sequence of times ti, • • where the are an index labeling the 
intervals. We can then specify one particular history by which interval was passed through at 
each time, labeling it by the sequence of indices ati, a 2 , ■ ■ ., which I will generally abbreviate 
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as a. Such a history includes all possible paths which pass through the given set of intervals 
at the given times. 

The decoherence functional is a functional on pairs of histories. The value of this func- 
tional on a pair of histories a and a' is given by 

D(a', a) = £ 5q> £ 5q 5(q' f - q f ) exp^(S[q'(t)} - S[q(t)])/h}p(q' , q ). (1.1) 

Here we are integrating over all paths which pass through the specified sequence of 
intervals. The functional S[q(t)] is the fundamental action. If ReD(a',a) = for a' ^ a, 
then the system is said to be decoherent, and obeys classical laws of probability. The diagonal 
elements D(a, a) are the probabilities of each history a. 

The simplest form of this type of history is that where the intervals are completely 
fine-grained in certain variables, and completely coarse-grained in others. We divide the 
coordinates q 13 into two groups, x 13 (henceforth known as the system coordinates), referred 
to collectively as x, and Q k (henceforth known as the reservoir coordinates), referred to 
collectively as Q. Our histories will then be complete trajectories x^it) for the system 
coordinates, while the reservoir coordinates will be neglected completely. 

It is then convenient to break the fundamental action of the system into several parts: 

S[q(t)} = S sys [x(t)} + S ies [Q(t)} - f tf V(x(t), Q(t)) dt, (1.2) 

J to 

where S' sys [a;(t)] is the action of the system, S res [Q(t)] is the action of the reservoir, and there 
is an interaction potential V(x, Q) between them. The decoherence functional is then 



4 



D[x'(t),x(t)} = exp^(S sys [x'(t)] - S sys [x(t)})/t} J ^Q' J SQ $(Q'(tf) - Q(t f )) 
x exp{i(S res [Q'(f)] - S ies [Q(t)\ 

- f tf (V(x'(t),Q'(t))-V(x(t),Q(t))) dt)/h\ 
xp(x' ,Q' ;x ,Q ). 



;i.3) 



II. LINEAR CASE 

The case of a system interacting linearly with a reservoir is a famous one, and has been 
treated by a number of people; quantum mechanically by Feynman and Vernon, and Caldeira 
and Leggett, classically by Zwanzig. For convenience, it is customary to make a number of 
simplifying assumptions: 

1. The reservoir variables Q k are harmonic oscillators, i.e., 

_ ft f jji . O 9 

S„[Q(t)] = £ / -(Q fe -co 2 k Q k2 ) dt. (2.1) 

2. The initial density matrix factors: 

p(x' , Q' ; x , Qo) = x(x' Q , x )(f)(Q' , Q ). (2.2) 

A similar assumption classically is to assume that the initial probability distribution of 
the reservoir coordinates is independent of the initial state of the system coordinates. 

3. The interaction V(x, Q) is bilinear: 

V{ Xl Q) = -Y,lkxQ k . (2.3) 

k 

I will generally assume that x is a single variable; multivariable systems are a trivial 
generalization, where the 7^ become matrices. 

We will relax these assumptions to a certain degree later on, but for now let us consider 
this case. The classical case is exactly solvable. In this, the equation of motion for the 
reservoir variable Q k is 
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^r(t) = -ulQ k (t) + ( lk /m)x(t). (2.4) 

This has a solution 

Q k (t) = Q k {t ) cos(u k (t - t )) + sm{u k {t - t )) + — f sm[u k {t - s)}x{s) ds. 

(2.5) 



We can then use this in the equation of motion for x: 



= (ttO (*) + F ^ + E — f sin M* - s)]x(s) ds, (2.6) 
V ox J i ~r mu k J to 



-0L\ ... yjLf 
^ mu k J to 

where F(t) is the Langevin force. In this case, it is clearly 



F(t) = E 7fc ( Qk(to) cos(cu fc (t - t )) + sm(u; k (t - tO))) . (2.7) 

If we assume that the Q have a thermal probability distribution initially, 



2vrA;T eXP V 2 

which, when averaged over an ensemble, gives 



= 0, (Q k2 ) = — 2 ,..., (2.9a) 
• 2 A-T 

(Q*> =0, (Q* > = —,... , (2.9b) 

7TZ 



then 



(F(0) = 0, (2.10a) 
W)^)> = E7, 2 (^)cos(^(t- S )). (2.10b) 

Let us compare this to the quantum results. It is interesting to first consider the system 
in isolation from the reservoir. In this case we would have 

D[x'(t),x(t)} = expfas^Wit)] - S sys [x(t)])/h)}xW ,x ). (2.11) 



If the action 5 , sys [z(t)] has the usual form 



S sya [x{t)\ = P L(x{t),x(t)) dt, (2.12) 



to 



then we can change variables 



X{t)= l -[x\t) + x{t% (2.13a) 
£(t)=x?(t)-x(t), (2.13b) 



and expand the phase in terms of £: 

s sys w(t)} - s sys [x(t)]= / ^mimm + f-nm^MO + cxo * 

J t oX 3X 

=jC'(-li (Jf(t) -* W)+ § wt) -* (t)) ) £W * 



%(X ,Xo)£(to) + 0(£ 3 ). (2-14) 
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So the Euler-Lagrange equation of motion appears in the phase of the decoherence functional! 

One should not put too much weight on this occurrence. This system is not decoherent; 
substantial interference can still occur between different possible trajectories. There is no 
particular reason to expect £ (t) to be small, so it is not correct to neglect higher-order terms. 
This system, on its own, is still essentially quantum-mechanical. It is not even quasiclassical. 

This still leaves the effects of the reservoir variables and interaction unaccounted for. Let 
us turn, then, to this portion of the decoherence functional. 



F[x'(t),x(t)} = I SQ' I 5Q5(Q'(t f ) - Q(t f ))exp^i(S Ies [Q'(t)] - S res [Q(t)} 

V(x'(t), Q\t)) - V(x{t), Q(t)) dt) /h} X {Qo, Qo) 
exp\iW[x'{t),x(t)]/h\. (2.15) 



to 



F[x'(t),x(t)] is termed the influence functional by Feynman and Vernon, and W[x'(t),x(t)} 
is the influence phase 0. In our simplified model, this is not difficult to evaluate exactly. 
It is generally assumed that the initial density matrix is in a thermal state. We quote the 
results of Feynman and Vernon: 



W[x'(t), x(t)} = - f tf dt I ds \x'{t) - x{t)} (k{t - s)x'{s) + k(t - s)x(s)) , (2.16) 

2 Jto Jt \ J 

where the real and imaginary parts of k(t — s) are 

k R (t-s)=J2 — MMt-s))i (2.17a) 
k rnu) k 

k^t - s) = J2 coth(huj k /kT) cos(uj k (t - s)). (2.17b) 
k mu k 



Changing to our variables X and £, we see that 

>/ 2 rt f rt 

k 



W[X(t),t(t)) = £ 77^ l f dtT ds{2Z(t)X(s) sin(c fc (t - s)) 
k lmu) k Jto Jt I 

coth cos(^(t - s))}. (2.18) 

Thus, we have a real term which is proportional to £(t) and an imaginary term which is 
proportional to £(i)£(s). The imaginary term is a double integral over a symmetric kernel 
whose eigenvalues are strictly non-negative; thus, for large £ the decoherence functional will 
be diminished by a decaying exponential 



exp 



dt / ds £(i)£(s) cos(a;fc(t — s)) 

to >/ to 



Since £ essentially measures how far you are from the diagonal of the decoherence functional, 
the off-diagonal terms tend to vanish and the system becomes decoherent. 

Furthermore, since large £ is suppressed, it now makes sense to discard terms of 0(£ 3 ). 
Thus we can now say 

S sys [x'(t)\ - S sys [x(t)} + W[x'(t),x(t)] = ~ f 1 dt [ tf ds £(t)ki(t - s)t(s) 

4 J to J to 

+ f tf dt£(t)<t)+0(e), (2.19) 

Jto 



where 

dt\dX J oX ^ muj k Jto 



If we compare this to Q2.6Q , we see that 



e(t) = 

is identical to the ensemble-averaged classical equation of motion. Note that the bath 

of harmonic oscillators acts as a retarded force on the system. In the limit as we go to 

a continuum of oscillator frequencies with a high cut-off, this retarded force becomes a 

dissipative term, i.e., a frictional force. In this limit, Caldeira and Leggett show that for a 

Debye distribution of oscillator frequencies, the influence phase becomes || 

rtt / . ikT \ 

W[X,$ = J (-2TXm + —Te(t))dt. (2.21) 

where T is the usual classical coefficient of friction, defined in terms of 7 and the cutoff 
frequency fl. See [|7|,|8| for details. 

We have seen that the real term of £] corresponds to the classical retarded or (in 

the limit) dissipative force. The imaginary term also has a classical analog. In the classical 
case, there is a random stochastic force F(t) given by ( |2.7| ), which ensemble-averages to zero 
(F(t)) = 0. As we see in (|2.10b|), however, the two-time correlation function of this force 



Im W[X(t),Z(t)} = f dt f ds (F(t)F(s))Z(t)Z{s). (2.22) 

Jtn Jtn 



does not vanish. As h —> 0, we get coth(huj/2kT) — > 2kT/hu. So the imaginary part of 
W[X, £] has the form 

r "dtf 

1 t J t 

Here we observe the subtle linkage between noise, dissipation, and decoherence. In 
interacting with the many degrees of freedom of the reservoir, the system loses energy. 
It also is subject to random jostlings from the reservoir oscillators. But one last, purely 
quantum-mechanical effect is that the state of the system is continually being "measured," 
and thus the various possible trajectories tend to decohere, at least on a scale large compared 
to h. Later we will see that even in situations where the classical noise vanishes, there is 
still quantum-mechanical noise. This arises essentially from the zero-point energy of the 
reservoir oscillators. 

We can straightforwardly generalize to the case where the potential is nonlinear in x, 
but still linear in Q. Suppose that 



V(x,Q) = -J2a k (x)Q k . 

k 

Here the influence phase is 

W[x'(t),x(t)} = £ t f dt fds {(a k (x'(t)) - a k (x(t))) 

k J to Jto I 

x(a k (x'(s))k k (t - s) - a k (x(s))k* k (t - s)) \, 



(2.23) 



(2.24) 



where 



k k (t - s) = 



sm(uj k (t — s)) +i coth(huj k /2kT) cos(u> k (t — s)) 



(2.25) 



2muJk 

We can again separate the real and imaginary parts, and change to variables X and £. We 
then get, to 0(£ 3 ), 

W[X(t),t(t)] = -7^E £ dt f to ds {a' k (X(t))a k (X(sMt)sm(Mt-s)) 

-i coth a' k (X(t))a' k (X(sMtMs) ™s(u; k (t - s))}. (2.26) 

Again, we see that the real term has the same form as the classical retarded force, which 
becomes dissipative in the limit of continuous frequencies and high cutoff. The imaginary 
term again corresponds to a double integral over the two-time correlation function of the 
classical stochastic force. It is strictly non-negative, and exponentially damps the decoher- 
ence functional for large £. 



III. NONLINEAR EXAMPLES 

The problem with potentials nonlinear in Q is that the path integral is no longer solvable 
in closed form. Thus, it is difficult to be certain that this correspondence with the classical 
equation of motion which holds in the linear case is truly universal. We can, however, 
consider weak couplings, and solve for the equation of motion using perturbation theory. 
We can then compare the classical perturbative equation to that derived from the influence 
functional. 
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A. Classical and Quantum perturbation theory 



Let us consider a system coupled to a bath of harmonic oscillators with a potential of 
the form 

V(x,Q) = -eJ2V k (x,Q k ), (3.1) 

k 

where Vk(x, Q k ) can be nonlinear in x and Q k . In general, such a problem cannot be solved 
exactly. However, if the coupling is weak (e << 1) then we can make a perturbation 
expansion, at least for reasonably well-behaved potentials. 
The total Lagrangian is 

U Ux, x, Q, Q) = L(x, x) + £ y ({Q k ) 2 - uj 2 k {Q k f^ - eV k (x, Q k ). (3.2) 

Let's suppose that the trajectory x(t) is known. Then the equation of motion for the kth 
harmonic oscillator is 

If we then write Q k as an expansion 



Q k (t) = Q k (t) + eQf(f) + e 2 Q k (t) + ... (3.4) 



and equate equal powers of e we get a series of equations 



^ = ~"lQl (3.5a) 



^ = + ^QJW^WO.OoW), (3.5c) 



etc., where we've Taylor-expanded V k (x, Qq + eQ\ + • • •) in powers of e. 
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Now we have equations for each Q\(t) in terms of the lower order functions. Notably, 
the lowest order equation is now a simple harmonic oscillator, and we can solve for it easily 
in terms of the initial conditions 

Qq{€) = A k cos(u k t) + B k sin(w fc t), (3.6) 

where A k = Q k \ to and B k = (l/uj k )(dQ k /dt)\ to . 

The higher-order equations are driven oscillators. We can solve for them exactly, match- 
ing initial conditions: 

1 f w hf ,. 



Qi(t) = / sinM* - s ))—±(x(s), Q k (s)) ds, (3.7a) 

muj k Jto oQ 

Q*(*) = — tsm(uj k (t-s))^(x(s),Q k (s))Q k (s) ds, (3.7b) 
mu) k Jt oQ z 

and so forth. 

Having found the motion of the harmonic oscillators in terms of x(t), we can now turn 
around and find the equation of motion for x. This is 

Q k {t) is the expansion that we solved for above, and it will depend on the earlier behavior 
of x, in general. Note that causality is strictly obeyed. This classical causality follows as a 
result of more fundamental quantum causality, as discussed by Gell-Mann and Hartle H]. 

We treat this same problem quantum-mechanically by trying to find the influence func- 
tional F[x'(t),x(t)} as a perturbation expansion. Assume that the reservoir starts in a 
definite initial state \a), with wave function <p a (Q). Then 

F a [x',x] = J 8{Q\t f ) - Q{t f )) exp ^(S ies [Q' (t)} - S TCS [Q(t)]) 
x 1 + ^ (v(x'(t),Q'(t)) -V(x(t),Q(t))j dt 

+ (^) 2 J'' j'(y(x\t),Q\t))-V(x(t),Q(t)^ 



'to Jto 

x ( V(x'(s), Q'{s)) - V(x(s), Q(s)) ) dsdt-t 



12 



MQ'(t m(Q(t )) SQ' 5Q 



l + j£(v aa (x'(t))-V m (x(t)f) dt 
~ (jpj jl{Va b {x\t))Vax'{s))e-^-^ - V ab {x{t))V ba {x\s))e-^-^ 

-y 6a (x / (i))y a6 (x( S ))e^(*- s ) + 4(i(t))^(3;( S ))e^Nj dsdt + --- 

i+ iiy L{xm{t)dt 

- [pi t,^ l\yL(x(t))v: b (x(sm(tns)cos(uj ba (t-s)) 



-2tV b ' a (X(t))V ab (X(sm(t) sin(uj ba (t - s)) ) ds dt + 0(e 3 ) + 0(£ 3 ). (3.9) 
Here we've defined the functions 

VU*) = (a\V(x,Q)\a) = I MQWa(QMx,Q) dr, (3.10a) 



V ba {x) = (b\V{x,Q)\a) = / MQWb(Q)V(x,Q) dr. (3.10b) 



In our case, we assume that the reservoir is a collection of harmonic oscillators initially 
in a thermal state. In this case, the states \a) become the ordinary Fock states \n) and the 
influence functional is 

F[x'(t),x(t)} = J2p nn F n [x'(t),x(t)}, (3.11) 

n 

where 

Pnn = ~ exp(-hu k /kT)^ exp(-n k hu k /kT). (3.12) 

B. Polynomial potentials 

We will specifically consider a potential of the form ( |3.1| ) where the individual potentials 
are polynomials in Q k . We will see that it is convenient to separate the even and odd terms: 

N N 

V k (x,Q k ) = ^a kl {x){Q k ) 2l+1 + J2b kl (x)(Q k f. (3.13) 

1=0 1=1 
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The ciki{x) and bki(x) are arbitrary functions of x, only assuming that the potential as a 
whole remains relatively well-behaved, integrable, etc. For convenience, I will drop the index 
k for the rest of this derivation. It should be understood that the final result is to be summed 
over all the oscillators, 

W[x'(t),x(t)} =J2W k [x'(t),x(t)}. (3.14) 

k 

From the equation ( |3.5aj ), we can write down the equations of motion for a classical 
oscillator Q(t) = Qo(t) + eQi(t) + . . .. We then plug in the solutions ( |3.6p and ( |3.7a| ) to get 



Qo(t) = A cos(ut) + B an(ut) = (3e lwt + (3*e~ luJ \ (3.15a) 



I "' I N N \ 

Q l {t) = — mn[ u {t-a)\ ^(2j + l)a i (a;( S ))Q^( S ) + ^ 2j'6 i (x( S ))Q^ 1 ( S ) ds, 

■"" \j () j | J 

(3.15b) 



etc., where (3 = (A — iB)/2. The equation of motion for x is then 



N N 



it (!) - (i) +«(E«jw«))or i w + E«i(«(«))#(o 

(N N 
£(2j + l)a}(x(0)g?(*) +E2j&;^W)Qo J " 1 W ] Qi(t) +0(e 3 ) 
3=0 i=l 



9x y 



+ e7i(t) + eV 2 (t) + 0(e d ) (3.16) 



We are interested in the ensemble-averaged equation. We can make use of the fact that 



(ran = s mn n\ [ ) . (3.17) 



So only the even terms contribute to the first-order component of the equation ( [3.161 ). Qq {t) 
is readily found then with a binomial expansion of 

(f3e iuJt + (3*e- iu>t f j = J2 f^l (J3e iut ) i (j3*e- iut f'~ i , (3.18) 



i=0 



/ 



n \ m 
A ) i\(n — i)V 
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yielding 



(3.19) 



The second order component is more complicated. Plugging expression (|3.15b| ) for Q\(t) 
into ( |3.16|) , doing a binomial expansion for the powers of Qo(t) and Qo(s), pairing e miu)t and 



-miujt 



terms, and ensemble-averaging gives us 



I rt r N N 

(t 2 (*)> = / EE sin [( 2 ^ + ~ s))C ijk (t, s) 

mu J t l k=Qi j=k 

N N 

-EE sin K 2 ^ - iM* - s)]c ijh (t, s) 

k=l i,j=k 
N N 

+ E E sm[2ku}(t - s)]D ijk (t,s) 

k=l i,j=k 

N N s 

-EE sin K 2k - 2)w(t - s)]D ijk (t, s) \ ds, 

k=l i,j=k 



where 



C ijk (t, s) = 4(x(t)) aj (x(s))(i + j)\{2i + l)(2j + 1) . 



(3.20) 



*i)(p-i){^T> (3 - 2ia) 



s) = 6((a;(t))6 i (a;(s))(i + j - l)!4ij 



'2z-l\ /2j - 1\ / feT V 4 "^ 1 



n « — k J \j — k J \ 2muj 2 J 
We can collect together and combine those terms with the same sine factor, to get 

1 rt( N N 

(r 2 (t)> = /EE sin t( 2fc + !M* " s)]^iifc(t, s) 

V V s 

+ E E s in [2^(t - s)]Fy fe (f, s) \ ds, 

k=0 i,j=k ' 



(3.21b) 



(3.22) 



where 

E ijk (t,s) = a' i (x(t))a j (x(s))(i + j)\ 



i+j 



2z + 1\ /2j + r 

(3.23a) 



i^(M) = 65(x(t))6,-(x( a ))(i + j - — 2 J 2fc(< + 



' 2. \ / 2j ^ 
i - k W - fc , 



(3.23b) 
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We are also interested in the correlation function (F(t)F(s)), where F(t) is the force due 
to the interaction with the reservoir. To second order this is 

( N N 

(F(t)F(s))=e 2 \Y, E 2cos[(2k+l)u(t- s)]G ijk (t,s) 

^k=0 i,j=k 
N N 

+ EE 2cos[2kto(t-s)]H ljk (t,s) 

k=l i,j=k 
N . 

+ E H ij0 (t,s)\ + O(e 3 ), (3.24) 

G^(t, a ) = aI(x(t))o;(x( 8 ))(< + j + l)!f^) +J+ ^k)^j-k)' (3 - 25a) 

H ijk (t,s) = b>(x(t))b>(x(s))(i +jV-[^J +3 (. !\) (/_^)- 0.25b) 

We can subtract off the average values to get 

{F(t),F{s)) = (F(t)F(s)} - (F(t))(F(a)), (3.26) 

where (F(t)) is the first order ensemble averaged force from ([3.19 ) . 

We can compare this result to that obtained from our quantum mechanical procedure. 
Suppose that the reservoir begins in a definite state \n). Then the influence functional is 
given by 

F n [X(t),m = 1 + ea nl [X(t),m + e 2 a n2 [X(t),^t)} + ■■■, (3.27) 
and in the thermal case by ( p. 11 ), 

F[X(t),£(t)] = 1 + «*i [*(*),£(*)] + e 2 a 2 [X(t), f (t)] + ■ ■ ■ = PnnF n [X{t)^{t)l (3-28) 

n 

where 

Oi[X(t),Z(t)] = J2PnnOCni[X(t),m\- (3-29) 
n 

The influence phase is then 
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W[X(t),£{t)] = -fohF[X(t),Z(t)] 
= -iheai[X(t),£(t)] 

-%ht 2 (a 2 [X(t),m ~ ~«i WU(Ol) +■■■■ (3-30) 

From ( |3.9| ), then, we see that we must find an expression for (m\r l \n). This will, in general, 
be a polynomial in n, for certain values of m, and zero for the rest. In comparing to the 
classical result, we need keep only the highest power of n, since the lower powers will be 
higher-order in fuo/kT as we let ft — > 0. This will be 

(m\r l \n) = (\ \(^—\ n l/2 + ■ ■ ■ , m = n + I - 2k, 2k < I, (3.31) 
V k I \ 2muj ] 



(m\r l \n) = (\ \(-^—\ m l/2 + • • • , m = n - I + 2k, 2k < I, (3.32) 
\kj \2muj J 

and zero otherwise. 

We can then use the fact that as h — ► 0, 

?"»-"' (3 ' 33) 



Thus, from equation ( |3.9|) we get 



a 1 {X(t),m=Y,Pnnl f f Y. h '^nt)){n\r 2n \n)m dt 

„ il J tn 



i ft f 2j\ ( kT 

3 

which agrees exactly with the first order term in the classical equation of motion ( p. 19 ) 
Similarly, we can calculate the second order term to get 

x , N N t t 



a 8 [X(fU(f)] = -=r{E E r I 2cos[(2k + l)u;(t-s)}G ijk (t,s)^(s)ds 

h l k=0i,j=k Jt ° J *o 



dt 



N N /•* 

+ E E / / 2cos[2fcu;(t- s)]H ijk (t, s)i{t)i(s) ds dt 

N N tf t 

+i E E / / sin[(2A; + l)u(t - s)]E ijk {t, s)£(t) ds dt 
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k= Hj =k J to J t 




i,j=l Jt " Jt ° 



(3.35) 



Here we've used the same definitions of E^k, etc., where the classical system variable x has 
become the quantum variable X. 

We can clearly see from this the exact correspondence with the classical equation of 
motion, at least to second order in e. The real part of H^[X(t), £(£)] is just an integral of 
the classical retarded force, just as in the linear case, and the imaginary part consists of a 
double integral 



note that the — (F(t)){F(s)) comes from subtracting al/2 from the second order term. 
Again, we note the non-negativity of this imaginary part; the presence of noise both makes 
the behavior unpredictable, and causes different trajectories to decohere. So we see that in 
perturbation theory, the nonlinear problem has exactly the same classical correspondence 
as the linear problem. 



Though the above discussion is fairly general, it leaves unexamined the far broader 
range of possible strong, nonlinear interactions, as well as the possibilities of non-oscillator 
reservoirs. This is, of course, a product of computational convenience, as it is very difficult 
to get analytical answers in other cases. Are there any arguments that can be made for 
more general systems? 

In any case where the action can be decomposed 




(3.36) 



IV. MORE GENERAL CASES 



S[x(t),Q(t)\ = S sys [x(t)\ + S ies [Q(t)\ + S int [x(t),Q(t)\, 



(4.1) 



it is possible formally to write the decoherence functional in the form 
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D[x(t),x'(t)] = exp -\ S sys [x(t)] - S sys [x'(t)} + W[x(t),x'(t)} 



(4.2) 



If we restrict ourselves, for the moment, to systems in a factorizable pure state, 



p(x',Q';x,Q) = *V)*(*)**(QWQ), 



then this influence phase is defined simply by (2.15) 



(4.3) 



exp {W[x'(t),x(t)]/h} = J 5Q'J 6Q5{Q'(t f )-Q(t f ))exp^S Tes [Q'(t)]-S ies [Q(t)} 

+S int [x'(t),Q'(t)] - S int [x(t),Q(t)]U*(Q' MQ ). (4.4) 



By bringing the integral over the final condition Qf, Q'j to the front, we can re- write this 
as a product of two path integrals: 



exp{iW[x'(t),x(t)]/h} = J I dQ f dQ' f 5(Q'(t f ) - Q(t f )) 

J 5Q'ex P ^S ICS [Q\t)] + S int [x'(t),Q\t)]^*(Q' ) 



x 



x 



SQexp -\ -S ies [Q(t)} ~ S int [x(t), Q(t)} \$(Q ) 



= J J 'dQ f dQ' f 5(Q'(t f ) - Q(t f ))$l, {t) (Q' f )<S> x{t) (Q f ) 

= J K'(t)(Qf)^(t)(Qf)dQf = I (4.5) 

where |3> x (t)) an d \&x'(t)) are the states that |<3>) will evolve into under the influence of the 
interaction, given the trajectories x(t) and x'(t), respectively. 

Clearly, (3> x '(t)|$x(t)) < 1> which implies equally clearly that lmW[x'(t),x(t)] > 0. So 
the non-negativity that we saw in the cases I and II above is generally true. This is also 
clearly the case for mixed states, since we can represent any mixed state as 



p(x',Q';x,Q) = J2Pi^')M^Wi(Q')HQ): 



(4.6) 



where 



(4.7) 
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so if the Fi[x'(t),x(t)] < 1, then clearly 



exp{tW[x'{t),x{t)]/h} = Y,&FiW(t),x{t)] < 1 



(4.8) 



and lmW[x'(t),x{t)] > still holds. Also, ImW[a;'(t),a;(f)] = for x'(t) = x{t). Thus, 
without assuming anything about the interaction or the reservoir, we see that there will be 
a maximum at = 0, and that the off-diagonal £(£) ^ terms will tend to be suppressed. 
This is not surprising, as one expects almost any sort of interaction with neglected degrees 
of freedom to result in the loss of phase coherence. However, it does show how these highly 
simplified models might actually demonstrate behavior important to the rise of classical 
physics from quantum mechanics in physical systems. 

For example, in considering quantum gravity, decoherence might arise from neglected 
gravitational degrees of freedom. The usual semi-classical treatment of quantum gravity, 
which omits the "back action" of mass-energy on the curvature of spacetime, cannot exhibit 
this effect. The weakness of the gravitational interaction would in general make it less 
important in causing decoherence than stronger forces, such as electromagnetism; but it 
might well become important in quantum cosmology. 

There are, of course, still questions. All that has been demonstrated is the non-negativity 
of lmW[x'(t),x(t)]. Can there not be zeroes for some choices of ^ 0? And how strongly, 
in general, are the off-diagonal terms suppressed? 

There can certainly be zeroes for nonzero in some cases. Indeed, if we consider the 
form of lmW[x'(t),x(t)} for the linear case 



(for a one-oscillator "reservoir" of frequency uj), there are an infinite number of choices of 
£(t) which make this zero. Thus, one cannot call this system truly decoherent. However, 
as the number of oscillator frequencies is increased, the number of possible choices of £ (t) is 
further and further restricted, so that as the reservoir becomes infinite only £(£) = remains. 
One would expect similar behavior in the more general case. While it is certainly possible 




t 



t 



(4.9) 



20 



to construct cases where lmW[x'(t),x(t)} has many zeroes even for a very large reservoir, in 
practice one expects lmW[x'(t), x(t)] > for x(t) ^ x'(t), as the degrees of freedom of the 
reservoir are increased. 

Similarly, the strength with which off-diagonal terms will be suppressed depends on the 
details of the system. However, one would expect that |$ x (t)) an d \&x'(t)) differ more in the 
case of strong interactions than small, and hence that (3> x '(t)|^x(t)) would be more strongly 
suppressed, in general. 

V. CONCLUSIONS 

It is clear that it is possible to define a "classical" equation of motion directly from the 
underlying quantum theory, and that, at least in many cases, this corresponds closely to the 
equation obtained from the classical theory. While correspondences of this sort have often 
been demonstrated in the past, never before has there been a rigorous, a priori technique 
for deriving them. 

Using the formalism of Gell-Mann and Hartle, we can now see classical physics as, very 
simply, a limit of the underlying quantum theory; and we can systematically determine, 
at least in principle, the deviations from strict classical equations due to quantum effects. 
Using the decoherence functional as a criterion for determining if an effect is experimentally 
observable, we can once and for all avoid the problem of collapsing the wave function; there 
is no longer any need for an independent "classical realm" of measurement. 
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After the completion of this research, I learned that Bei Lok Hu, Juan Pablo Paz, and 
Yuhong Zhang had studied a very similar class of nonlinear brownian systems more or less 
simultaneously with me ||10|| . While their study is from a considerably different point of 
view, with very different goals, being chiefly concerned with deriving master equations and 
treating the thermodynamics of these generalized systems, their results overlap mine to a 
certain extent. 
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